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Abstract 

The apparent dependence of nuclear matter observables on off-shell prop- 
erties of the two-nucleon potential is re-examined in the context of effective 
field theory (EFT). Finite density (thermodynamic) observables are invari- 
ant under field redefinitions, which extends the well-known theorem about 
the invariance of S-matrix elements. Simple examples demonstrate how field 
redefinitions can shift contributions between purely off-shell two-body inter- 
actions and many-body forces, leaving both scattering and finite-density ob- 
servables unchanged. If only the transformed two-body potentials are kept, 
however, the nuclear matter binding curves will depend on the off-shell part 
(generating "Coester bands"). The correspondence between field redefinitions 
and unitary transformations, which have traditionally been used to generate 
"phase-equivalent" nucleon-nucleon potentials, is also demonstrated. 
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I. INTRODUCTION 



Over thirty years ago, it was observed that "phase-equivalent" nucleon-nucleon potentials 
gave different predictions when used to calculate the energy per particle of nuclear matter 
£j]-§]|- Phase equivalent means that the potentials predict identical two-nucleon observables 
(phase shifts and binding energies); families of such potentials can be constructed via unitary 
transformations ||. The finite-density discrepancies were attributed to the differing off- 
shell behavior of the nucleon-nucleon (NN) amplitude or T-matrix, and were related semi- 
quant it at ively to the distortion of the two-body relative wave function in the medium (the 
"wound integral") |!]]. It was proposed that comparisons of nuclear matter calculations could 
help to determine the correct off-shell dependence [f|.|6|]. 

From the modern perspective of effective field theory (EFT), these results demand a 
different interpretation, since off-shell effects are not observable. Phase-equivalent EFT 
potentials should give the same result for all observables (up to expected truncation errors) 
if consistent power counting is used; differences imply that something is missing. In this 
paper, we show how this issue is resolved by considering field redefinitions at finite density 
in an effective field theory. 

There is an extensive literature, primarily from the sixties and seventies, on the role 
of off-shell physics in nuclear phenomena (see, e.g., Ref. || and references therein). This 
includes not only few-body systems (e.g., the triton) and nuclear matter, but interactions 
of two-body systems with external probes, such as nucleon-nucleon bremstrahlung and the 
electromagnetic form factors of the deuteron. The implicit premise, which persists today, is 
that there is a true underlying potential governing the nucleon-nucleon force, so that its off- 
shell properties can be determined. Indeed, the nuclear many-body problem has traditionally 
been posed as finding approximate solutions to the many-particle Schrodinger equation, 
given a fundamental two-body interaction that reproduces two-nucleon observables. 

An alternative approach to nuclear phenomena is that of effective field theory. The EFT 
approach exploits the separation of scales in physical systems l?j-|T3f. Only low-energy (or 
long-range) degrees of freedom are included explicitly, with the rest parametrized in terms of 
the most general local (contact) interactions. Using renormalization, the influence of high- 
energy states on low-energy observables is captured in a small number of constants. Thus, the 
EFT describes universal low-energy physics independent of detailed assumptions about the 
high-energy dynamics. The application of EFT methods to many-body problems promises a 
consistent organization of many-body corrections, with reliable error estimates, and insight 
into the analytic structure of observables. EFT provides a model independent description 
of finite density observables in terms of parameters that can be fixed from scattering in the 
vacuum. 

A fundamental theorem of quantum field theory states that physical observables (or more 
precisely, S-matrix elements) are independent of the choice of interpolating fields that ap- 
pear in a Lagrangian Jl4],[l5| . Equivalently, observables are invariant under a change of field 
variables in a field theory Lagrangian (or Hamiltonian) . This "equivalence theorem" holds 



for renormalized field theories, although there are subtleties with some regularizations [TL6 



In an EFT, one exploits the invariance under field redefinitions to eliminate redundant terms 
in the effective Lagrangian and to choose the most convenient or efficient form for practical 



calculations [17-22]. Since off-shell Green's functions and the corresponding off-shell am- 



plitudes do change under field redefinitions, one must conclude that off-shell properties are 
unobservable. 

Several recent works have emphasized from a field theory point of view the impossi- 
bility of observing off-shell effects. In Refs. [f^UM], model calculations are used to illus- 



trate how apparent determinations of the two-nucleon off-shell T-matrix in nucleon-nucleon 
bremstrahlung are illusory, since field redefinitions shift contributions between off-shell con- 
tributions and contact interactions (see Sec. [TIT]). Similarly, it was shown in Ref. that 
Compton scattering on a pion cannot be used to extract information on the off-shell behavior 
of the pion form factor. The authors of Refs. pE| , |27|] emphasized the nonuniqueness of chiral 
Lagrangians for three-nucleon forces and pion production. Field redefinitions lead to differ- 
ent off-shell forms that yield the same observables within a consistent power counting. In 
Ref. [2j|, an interaction proportional to the equation of motion is shown to have no observ- 
able consequence for the deuteron electromagnetic form factor, even though it contributes 
to the off-shell T-matrix. 

In systems with more than two nucleons, one can trade off-shell, two-body interactions 
for many-body forces. This explains how two-body interactions related by unitary trans- 
formations can predict different binding energies for the triton [^] if many-body forces are 
not consistently included. These issues are discussed from the viewpoint of unitary trans- 
formations in Refs. [3(J and [31]. An EFT analysis of the triton with short-range forces was 
performed by Bedaque and collaborators. They showed that a three-body counterterm is 
needed in the triton for consistency JHJ . Varying this counterterm generates a one-parameter 



family of observables, which explains the Phillips line. 

Here we consider fermion systems in the thermodynamic limit. We extend the equiva- 
lence theorem to finite density observables, or thermodynamic observables in general, and 
show the consequence of field redefinitions. Since the thermodynamic observables are con- 



veniently expressed in a functional integral formalism, we adapt the discussion in Ref. pi 
which provides a particularly transparent (if not entirely rigorous) account of the equiva- 
lence theorem for S-matrix elements and how it can be applied in EFT's, based on functional 
integrals. We show how field redefinitions applied to an EFT Lagrangian with only on-shell 
vertices generate new off-shell vertices, which by themselves change predictions for an in- 
finite system. However, many-body vertices are also induced, and if included consistently, 
will precisely cancel the off-shell parts. In the end, all calculations agree, with no off-shell 
ambiguities. 



In section II, we review the argument in Ref. ||21||, extend it to thermodynamic observ- 



ables, and show how field redefinitions can reduce an effective Lagrangian to a canonical 
form. In section III, we illustrate the effects of field redefinitions with a simple example 
applied to the dilute Fermi gas of Ref. | 3~3fl . The connection to traditional treatments is 
elucidated in section IV. Finally, we summarize in Section V. 



II. FIELD REDEFINITIONS 



The equivalence theorem is typically stated as the invariance of the S-matrix for asymp- 
totic scattering states under a change of field variables. Intuitive diagrammatic proofs, as 
in Refs. |34|,[35|], are based on the reduction formula for calculating S-matrix elements from 



on-shell Green's functions. The proofs show that field redefinitions change only the wave 
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function renormalization, which is divided out in the reduction formula. The extension to 
observables in the thermodynamic limit, such as the energy density of an infinite system of 
fermions at zero temperature, is not obvious since there are no asymptotic single-particle 
states. Here we make that extension for nonrelativistic systems and discuss the conse- 
quences for effective field theories. The generalization to covariant effective field theories is 
st raight forward . 

A rigorous proof of the equivalence theorem must include a careful discussion of renor- 
malized Green's functions and amplitudes, and the impact of different regulators. Such 
a proof, based on the renormalized quantum action principles, is given in Ref. jT5[ and 



extended to a heavy-fermion expansion (in the context of heavy quark effective field the- 
ory) in Ref. ||19|| . The new issues in considering thermodynamic observables, however, do 
not change the subtle technical issues in the proof, so we build instead on a less rigorous 



but more intuitive path integral discussion of the equivalence theorem |2T| to illustrate the 
broad concepts and consequences. The path-integral formalism is convenient, since one can 
generate Green's functions and subsequently S-matrix elements, as well as thermodynamic 
observables, in the same framework. 

We will only consider heavy (nonrelativistic) boson or fermion fields explicitly; the gen- 
eralization to include long-range bosonic fields (e.g., such as the pion) is straightforward. 
The procedure that leads from an underlying relativistic theory to an effective Lagrangian £ 



in a heavy-particle formulation is outlined in Ref. |K|. We work in the rest frame, so there 
are no non-trivial four-velocities associated with the field variables. In general, £ contains 
all terms consistent with the usual space-time symmetries as well as any symmetries of the 
underlying theory. These include, in general, higher-order time derivatives, which can be 
eliminated as described below. The path integral form of the generating functional is written 
schematically as: 

Z[^,i] = /p^t^ei/^W^? ; (!) 



where ip and ijy are Grassmann variables for fermions, £t and £ are external Grassmann 
sources. By taking functional derivatives of Z with respect to the sources, we generate the 
n-body Green's functions and then S-matrix elements via a reduction formula. 

Here we follow the discussion of the equivalence theorem in Ref. ||21|| . The basic obser- 
vation is that if we change field variables in the functional integral, predictions for given 
external sources should be unchanged. Consider a transformation of ip of the form: 

ij;{x) — >^{x) +r]T[ijj,^} , (2) 

where T[ip,ip'] is a local polynomial in ip{x), ip'(x), and their derivatives, and r] is an 
arbitrary counting parameter (to be used below). The transformation given by Eq. (^|) 
is local and respects the symmetries of £. There are no space-time dependent external 
parameters, i.e. the x-dependence of T[ip, ip'] comes from the fields alone. Specific examples 



of such transformations are given in Eqs. (18) and ([50 



Within the path integral, three types of change result from Eq. (0): 

1. the Lagrangian changes: £ — > £'; 

2. there is a Jacobian from the transformation; 
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3. there are additional terms coupled to the external sources £ and £t. 



The statement of the equivalence theorem is that changes 2 and 3 by themselves do not 
affect observables, so that one is free to change variables in the Lagrangian alone. By 
introducing ghost fields, the Jacobian from 2 can be exponentiated. As detailed in Ref. plj . 



new terms involving ghosts either decouple or contribute pure power divergences in ghost 
loops. The latter are either zero in dimensional regularization or can be absorbed without 
physical consequence into couplings in C Therefore the Jacobian has no observable effect 
The extra terms in 3 do change the Green's functions but not the S-matrix. The proof is 
analogous to the usual diagrammatic proof ||21|| . 

To extend the discussion to finite density and/or temperature, we consider a Euclidean 
path integral formulation of the partition function Z and thermodynamic potential Q: 

Z = exp(-/3fl) = Tr{exp[-/3(F - /iN)]} , (3) 

where (3 is the inverse temperature, /i is the chemical potential, H is the Hamiltonian, 
and N is the number operator (e.g., the baryon number for nuclear systems). The hat 



indicates second quantized operators that act in Fock space. We direct the reader to Ref. |37 
for details on the derivation and interpretation of nonrelativistic path integrals at finite 
temperature and chemical potential. The connection between H and £ is complicated, but 
we do not need to actually derive the connection, since we assume a general form for C. 
However, we need to be more explicit with N, which is not in general equal to the spatial 
integral of ijrijj. Instead, the appropriate operator should be derived using Noether's theorem 
applied to C. 

The integrations in the path integral corresponding to Eq. ([3]) will be over complex 
variables satisfying periodic boundary conditions in imaginary time for bosons and over 
Grassmann variables satisfying antiperiodic boundary conditions for fermions. If there is 
no chemical potential, the extension of the equivalence theorem to finite temperature is 
immediate, because the finite temperature boundary conditions are trivially preserved by 
field redefinitions of the form in Eq. @) and the remainder of the argument is unchanged. 
Thus the thermodynamic quantities derived from the partition function (e.g., free energies 
and specific heats) will be invariant under field redefinitions for all (3. Since there are no 
external sources, we are free to introduce linearly coupled sources to construct a perturbative 
expansion after the field redefinitions. 

A nonzero chemical potential will be multiplied by new terms after the field redefinition. 
The conserved current is associated with a global U(l) symmetry transformation 

V>(z) — ► e-^iff(x) and ^{x) — ► e l ^{x) , (4) 

under which £ is invariant. We can conveniently identify the Noether current by promoting 
4> to a function of x and considering infinitesimal transformations with f36| 



1 Note that the field redefinitions Eq. (|2|) are different from the transformations used to explore 
anomalies, where contributions from the Jacobian are critical [36|. The latter transformations 
involve space-time dependent external parameters. 
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£ — • (5) 
Then the charge density p is given by 

p = ^; <j>(x)] = Fty, ^] ■ (6) 

The possibility of higher-order time derivatives or time derivatives in non-quadratic terms 
in £ means that, in general, F[ip,ip^] ^ ip^ip. 
Under the transformation Eq. (0), 



C[if>, V f ] — > £'[V>, = + 77T, V f + 

Sip 



while 



ip 1 } — > F'[ip, = F[if> + r]T, ^ + r]T f ] 



= F+(r / T- + H.c.)+... , (8) 

where we have written the first terms in the functional Taylor expansions of Eqs. (|7], |8|) and 
H.c. denotes the Hermitian conjugate. If T contains no time derivatives (as in the examples 



in Sections and [TV]), then using the fact that 8/8(d t (j)) commutes with 5 /Sip, we see 



order-by-order in 77 that 

p = iJr^Z'W'^W = F '^M ■ ( 9 ) 

5(d t <p) 

So the term multiplying p in the transformed path integral is the new charge density. If T 
contains time derivatives, p' will contain additional terms proportional to the equations of 



motion, which give zero contribution [19,21 



At zero temperature, we can work without an explicit chemical potential, but instead 
specify the density through boundary conditions in the noninteracting propagator (see 
Sect. |H). Using the new charge density operator, we find the density is unchanged by 
the transformation. This is illustrated explicitly for a special case in the Appendix. 

For an effective field theory, invariance under nonlinear field redefinitions has practical 
and philosophical implications. On the practical side, because the EFT already contains all 
terms consistent with the underlying symmetries (which are preserved by the field variable 
changes we consider), there is no penalty for making such transformations: they just change 
the coefficients of existing terms. Therefore, we can use this freedom to remove redundant 
terms and define a canonical form of the EFT Lagrangian. This can simplify the calculation 
or improve its convergence. One possibility is what Georgi calls "on-shell effective field 
theory," which features only those vertices that are non-zero for on-shell kinematics (of free 
particles) [|1^]. This was the choice in Ref. [33[|, which treated dilute Fermi systems. 



An effective Lagrangian is organized according to a hierarchy, which we label by 77 pl |: 

00 

£ = Y,V n £n- (10) 



ra=0 
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For example, we have rj oc 1/A for a natural field theory as in Ref. [Q. Define Co to be: 



idt + 



2M 



As a result of the transformation Eq. (|j), £ — > with 

^1 



£' = £0+ (r]T[ip,^) 



id t + 



2M 



1> + H.c. + £ /f 4 



(11) 



(12) 



n=l 



Thus, at leading order in 77, we generate a term proportional to the free equation of motion. 
The corresponding interaction term will be proportional to the inverse free propagator, and 
thus will vanish on shell. So we can change off-shell Green's functions (e.g., the two-nucleon 
off-shell T-matrix) by field transformations without changing observables. The consequence 
is that one can change the appearance of C. One may find that better convergence of an EFT 
expansion is correlated with a particular off-shell dependence, but the idea of measuring 
the "correct" off-shell behavior has to be dismissed. These implications are extended to 
conventional approaches to nuclear phenomena by considering them as implementations, 
albeit incomplete, of the EFT approach. The discussion can be directly extended to a finite 
system with a single particle potential included in Co (and in the noninteracting propagator). 

We are also able to remove redundant terms, which was implicit in the choice of the 
"general" effective Lagrangian in Ref. p3f. Time derivatives other than the quadratic term 



can be systematically eliminated by field redefinitions of the form in Eq. (fj) by exploiting 
the EFT hierarchy. In particular, the lower-order (in rf) terms can be removed at the cost 
of changing the coefficients of higher-order terms in 77. For example, suppose we have the 
Lagrangian 



id t + 



2M 



^ + r ll ^(td t ) 2 i; + r ]2 ^^( l d t )i;+ H.c.) + 0( V 2 ) , 



(13) 



where V The field redefinition 



2 



2M 



^ — ^ + 1 W f ) + (I - 4 ^)^ f » 



results in the new Lagrangian 



id t + 



2M 



(14) 



(15) 



with no time derivatives outside of Co to order 77. Once we remove terms up to some order 
in rj, further transformations will not reintroduce them because they involve the next higher 
order .0 



'For another example of such a field redefinition see Ref. [38| 
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Equation flT2|) shows that different choices of effective Lagrangians, related by field redef- 
initions, can differ by off-shell contributions. The equivalence theorem guarantees that these 
must ultimately be canceled in observables (or vanish individually). In a sense, that makes 
any explicit demonstration "trivial" . The realization of the equivalence theorem, however, is 
often counterintuitive. Approximations in a truncated theory have to be consistent with this 
realization. A common source of confusion arises in discussions that assume only two-body 
vertices are kept in many-body calculations. Therefore, we turn to some concrete examples. 



III. PERTURB ATIVE EXAMPLES 
A. Transformation of the Lagrangian 

In this section, we illustrate the connection between off-shell two-body terms and three- 
body terms using a simple field redefinition. In particular, the effect of an off-shell vertex 
generated by a field redefinition is exactly canceled by a three-body term generated by the 
same transformation. We demonstrate this cancellation explicitly for the process of 3 — >• 3 
scattering in the vacuum and for the energy density of a dilute Fermi system at finite density. 
For simplicity, we consider a natural EFT for heavy nonrelativistic fermions of mass M that 
has spin independent s-wave interactions whose strength is correctly estimated by naive 
dimensional analysis. Such a theory is discussed in detail, e.g., in Refs. 

For a natural theory all effective range parameters are of order 1/A, where A is the 
breakdown scale of the theory. For momenta k <C A, all interactions appear short-ranged 
and can be modeled by contact terms. We consider a local Lagrangian for a nonrelativistic 
fermion field that is invariant under Galilean, parity, and time-reversal transformations: 



C = ip 



td t + 



— yi 

V V - + ^[(V'VO^V 2 ^) + H.c.l + . . . , (16) 



2M 



16 



where V is the Galilean invariant derivative and H.c. denotes the Hermitian 

conjugate. Higher-order time derivatives are not included, since it is most convenient for our 
purposes to eliminate them in favor of gradients, as described in Sec. |I| Naive dimensional 
analysis estimates Cn to be roughly 4ir/(MA 2t+1 ), which implies a perturbative expansion 
in k/A. Co and C2 can be related to the s-wave scattering length a and effective range r e 
via a matching calculation. This leads to the relations 

Co = ^ and C 2 = C ^. (17) 

We stress that we choose this simple example to make the cancellation of off-shell ef- 
fects in observables as transparent as possible. As follows from the general discussion in 
the previous section, analogous cancellations occur in more complicated theories that in- 
clude long-range physics, relativistic effects, spin-dependent interactions, nonperturbative 
resummations and so on. 

We generate a new Lagrangian £' by performing the field transformation 
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Pi, 



K,Pi) 





+ 



1^2, P 2 j 



K,p' 2 ) 




+ 




+ 



FIG. 1. Two-body vertices generated by £'. The corresponding Feynman rules are given in 
Eq. pi). 





+ 




+ 



IV 



(v) 



FIG. 2. Three-body vertices generated by £'. The corresponding Feynman rules are given in 
Eq. (P). 



on L. The factor 1/A 3 is introduced to keep the arbitrary parameter a dimensionless. The 
additional factor of Air is introduced for convenience. With this choice, the induced terms 
have the size expected from naive dimensional analysis if a is of order unity. We obtain for 



£ = C 



Ana 



+ 



+ 



A 3 
AnaC 2 

Ana 



2(^V)(^) f (^V 2 ^) 



A 3 



+ H.c. 



+ 4- and higher-body terms + 0(a' 



(19) 



where 4- and higher-body terms as well as terms of 0(a 2 ) have been omitted. 

The two-body vertices generated by the transformed Lagrangian C are shown in Fig. |l|. 
The vertices (i) and (ii) were already present in the untransformed Lagrangian C, while 
(iii) is a new energy-dependent off-shell vertex of 0(a), which vanishes identically if all 
lines are on-shell. However, £' also contains two new three-body vertices of 0(aCo) and 
O(aC2).0 These vertices are illustrated in Fig. |[ The spin projection of a particle with 
energy Wj(u^) and momentum Pi(p'j) is denoted by s^s'j); we abbreviate <5 S . S ^ = 5ij>. Note 
that Ui p\j (2M) unless the particle is on shell. The Feynman rules for all vertices are 
listed below. 



These vertices are already present in a general version of £, in which case the new contributions 
change the corresponding coefficients. The cancellation discussed here is unchanged. 
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Two-body vertices: 



(i) -iC S 2 

.C; 



ii) _i y{(Pi -P2) 2 + (p'i -P 2 ) 2 }^ 
4ira ( 1 



(iii) 



A 3 



UJ 1 +U 2 + oj x 



UJr, 



2M 



P? + P2 + p'i + pI 



S 2 



(20) 



Three-body vertices: 



iv 



(v) 



A 3 
Area C 2 



12C 5 3 



Pi - P2j 



Perm{l',2',3'} 



+ (Pl - PS)' + (Pi - Pi)' + (P2 - PS)' + (P2 - P' 2 ) 2 



J ^2 



+ (P3 " P' 3 ) 2 + (P'i " P 2 ) 2 + (P'l - P' 3 ) 2 + (P2 " P' 3 



(21) 



The sum in (v) runs over all permutations of {1', 2', 3'} and is the totally antisymmetric 
Levi-Civita tensor with eyyy = —tyvz 1 = 1- Note that vertex (v) can be rewritten in terms 
of the total momentum P = pi + p 2 + P3 as 

Ana„ ^ „ , „ r 2 

3P - Pl • P 1 - P2 ■ P 2 - P3 • P 3 



[V) -i—r^-C 2 Sty 622^331 €i'2'W 

Perm{l',2',3'} 



(22) 



The spin functions 5*2 and 5*3 are given by 

5 2 — ^ll"^22' — ^12'^21' > 

53 = <5ll'(#22'£33' — ^23'^32'J + $L2' (^23'^31' — ^21"^33') + ^13' (^21'^32' — <^22'^31') 



(23) 



for spin-1/2 fermions and by S 2 = 2 and S3 = 6 for spin-0 bosons. The propagator for a 
particle line with energy u and momentum q is simply 



(24) 



lu - q 2 /(2M) + ie 

As discussed in Sec. [0], the Lagrangian £' has to give the same results for observables as 
the original Lagrangian C. For this to happen, certain combinations of diagrams involving 
vertices (iii), (iv), and (v) have to cancel. As a, C , and C 2 are independent parameters, this 
has to occur at every order in a, Co, and C 2 independently. In the following two subsections, 
we illustrate this cancellation for 3 — ► 3 scattering in the vacuum and for the energy density 
at finite density. 



B. Cancellations at Zero Density 

As mentioned above, all diagrams within one order of a have to cancel identically. Fur- 
thermore, the power counting of the original, untransformed theory carries over. The sim- 
plest example is the tree-level 2^2 diagram with vertex (iii) which vanishes identically 
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(a) (b) ( C ) 

FIG. 3. Additional diagrams for 3 — > 3 scattering generated by £ at (D(aCb). The external 
lines are labeled as in Fig. ^. 

when all external legs are on shell so that the cancellation at 0(a) is trivial. Similarly, it 
is straightforward to show that all diagrams at 0(aCQC l 2 ) that are bubble chains including 
k vertices of type (i) and I vertices of type (ii) plus one vertex of type (iii) vanish as well. 
This is most easily seen in dimensional regular izat ion, where all vertices are on shell. 

In the following, we consider 3^3 scattering, where cancellations between different 
diagrams occur. We concentrate on the case of spin- 1/2 fermions. The case of spin-0 bosons 
is analogous but somewhat simpler, because the spin function is trivial. At O(aCo), we have 
the three tree-diagrams shown in Fig. ^. The external lines are labeled as in Fig. ^. Using 
the Feynman rules from Eqs. (p0[)-(2^), we find for diagram 0(a), 



-M 3 ( a ) = i-fcfCo (5n'5 2 2'^33' - ^ll'^23'^32' ~ #12' #21' #33' + S 13 > 6 2 l> S 32 ') , (25) 

where the internal propagator was canceled against the corresponding dependence of the off- 
shell vertex. At this point, we have to take into account that the cyclic permutations of the 
incoming and outgoing lines in diagram [5](a) generate topologically distinct diagrams that 
have to be included as well. Consequently, diagram 0(a) really represents nine diagrams, 
which when taken together generate the full spin function 5*3. After summing over the cyclic 
permutations of the incoming and outgoing lines, we obtain 

4:7TQ 

-M 3 (a) = Q^CoSs = Mg h) . (26) 

Evaluating diagram 0(b) gives exactly the same result as diagram @(a). Finally, diagram 
0(c) is simply given by vertex (iv). Adding the contribution of all diagrams at 0{olCq) we 
find 

■M 3 ( a ) + A%b) + Mjfa = , (27) 

and to 0(olCq) the off-shell dependence introduced by vertex (iii) cancels identically. A 
similar cancellation occurs for the diagrams at OiptC'i) where the Cq vertex in Fig. is 
replaced by a C 2 vertex and vertex (iv) is replaced by vertex (v) . The corresponding diagrams 
are shown in Fig. [|. Due to the momentum dependence of the vertices, the calculation is 
somewhat more tedious but straightforward. After summing over the cyclic permutations 
of the external lines, we find 

A*4 W + -M| (b ) = -Mg c) , (28) 

where A1|( c ) is simply given by vertex (v) in Eq. (|2T|). As a consequence, all diagrams at 
first order in a cancel identically and the S-matrix elements for 2 — > 2 and 3 — > 3 scattering 
are unchanged by the field redefinition Eq. fll8f) . 
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-^C^ + ^^^^ + 




(a) (b) (c) 

FIG. 4. Additional diagrams for 3 — > 3 scattering generated by £' at 0(062). The external 
lines are labeled as in Fig. ^. 

C. Cancellations at Finite Density 

We now proceed to finite density and consider a dilute fermion system at zero temper- 
ature. We again concentrate on the case of spin-1/2 fermions. The extension to include 
higher spins or isospin is straightforward and would only change the spin degeneracy factor 
g. It is not necessary to explicitly introduce a chemical potential in this case. Note that 
the number operator changes with the transformation. The invariance of the density for a 
uniform system is illustrated in the Appendix. 

The natural EFT at finite density is discussed in detail in Ref. |J3]. Here, we only repeat 
the Feynman rules for the energy density: 

1. Assign nonrelativistic four-momenta (frequency and three-momentum) to all lines and 
enforce four-momentum conservation at each vertex. 

2. For each vertex, include the corresponding expression from Eqs. (^, [H]). For spin- 
independent interactions, the two-body vertices have the structure {5\^5ps + 8\sfip~f), 
where X, (3 are the spin indices of the incoming lines and 7,<5 are the spin indices of 
the outgoing lines.f] For each internal line include a factor iGo(fc)A 7 , where k = (k , k) 
is the four-momentum assigned to the line, A and 7 are spin indices, and 

~ U-S^U + U ( v(2mU ) • m 

3. Perform the spin summations in the diagram. In every closed fermion loop, substitute 
a spin degeneracy factor — g for each 8\\. 

4. Integrate over all independent momenta with a factor / d A k/ (27r) 4 where d A k = dko d 3 k. 
If the spatial integrals are divergent, they are defined in D spatial dimensions and 
renormalized using minimal subtraction as discussed in Ref. . For lines originating 
and ending at the same vertex (so called tadpoles), multiply by exp(ik r]) and take 
the limit 77 — > + after the contour integrals have been carried out. This procedure 
automatically takes into account that such lines must be hole lines. 



1 



For diagrams with no external legs it is more convenient to take care of the antisymmetrization 
in the spin summations. This is discussed in rule |3|. 



12 




FIG. 5. Finite Density: contribution to the energy density at 0(a). 



+ 0«£0 + 





(a) (b) (c) 

FIG. 6. Finite Density: contributions to the energy density at 0{olCq). 

5. Multiply by a symmetry factor ij (5'n|=2 X (^) m ) where S is the number of vertex per- 
mutations that transform the diagram into itself, and m is the number of equivalent 
/-tuples of lines. Equivalent lines are lines that begin and end at the same vertices 
with the same direction of arrows. 

When an off-shell vertex cancels a propagator, a new tadpole can be generated or a tadpole 
loop can be shrunk to a point. We need an additional rule to deal with this situation. 

6. New tadpoles require a convergence factor as described in rule 4. Tadpole loops shrunk 
to a point are zero. 

At 0(a) there is only one diagram, which is shown in Fig. |5|. Using the Feynman rules 
from above and the previous section, we find 

■ 4 |?{2*o + 2*-^(2k 2 + 2 q *)}. (30) 



XJ- 



After performing the contour integrals, go an d ko pick up the poles from the corresponding 
propagators, and we find £5 = 0. Thus the diagram at 0(a) vanishes identically, as in the 
vacuum. 

We now proceed to 0((xCq), where we have the three diagrams shown in Fig. |6]. Using 
the Feynman rules, we obtain for diagram |6|(a) 

i Attcx f d^p f d^k f d^q — ~ 

£ 6{a) =2<7(<7-l)— J j—y A J j—^ J ^-^G (p)G (k)G (p-q)G (k + q) 

x {2£;o + 2p - (p 2 + k 2 + (k + q) 2 + (p - qf 



Ana. 

2g(g - l)i-^C7 



(2tt) 4 UV ; 



2 g(g - 1)^C _^I_ (31) 

yvy ) A a °2167r 6 ' v ; 



where the off-shell vertex was canceled against the propagators with appropriate momenta. 
Interestingly, the diagram ||(b) does not vanish separately even at T = 0, as one would 
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+ 




+ 




(a) (6) (c) 

FIG. 7. Finite Density: contributions to the energy density at OiaC'i). 



naively expect. The reason is that the off-shell vertex cancels one of the propagators in 
the middle loop with two propagators. Consequently, the product 9(p — hp) 9 (hp — p) corre- 
sponding to a particle-hole pair with equal momentum, which usually makes such diagrams 
vanish at T = 0, does not arise. In particular, we have 



6(b) 



-9(9 



^ 2 i Ana 



d 4 k 



d 4 q f d A p 



x \ 2k + 2p - 



2M 



A 3 / (27r) 4 7 (2vr) 4 y (2vr) 



2p 2 + 2k' 



-2*7(0-1) 



2 47T« k F 



A 3 u 216tt 6 ' 
Finally, we obtain for diagram |6|(c) 



^ k °e i ^G {k)G (q) (G (p)Y 



(32) 



1 Attc\ 

£ 6 (c) = -9(9 - 1)(9 - 2)— (-i)12C„^ 



d A k 



iG {k) 



A 3 u 216tt 6 



(33) 



Adding all contributions to the energy density at O(aC ), the sum of £^( a ) and £|(b) builds 
up exactly the same spin degeneracy factor as occurs in £|( c ), and we have 

S$(a) + %b) + %c) = . (34) 

The corresponding diagrams at 0{aC<2) are shown in Fig. |7|. After their evaluation, we find 

^7(a) = 9{9 ~ lJ^^Ca 



A 3 ' 180vr 6 ' 
1 , 1 o Aira „ kh 1 



C 2 



1807T 6 



47ra 



^ (c) = -^-l)(^-3)— -C 2 



fc 
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A 3 z 180tt 6 



(35) 



where the sum of £|( a ) and £^b) again builds up the same spin degeneracy factor as in 8^ 
and the net contribution of all diagrams at 0{aC<i) is zero. In contrast to the contribution 
at O(aC ), the induced three-body term with derivatives generates a factor of g — 3 instead 
of g — 2 and does not vanish separately for spin- 1/2 fermions. In conclusion, the contribution 
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FIG. 8. Energy density for a dilute Fermi system with spin degeneracy g = 4 through 0(k F ) 
as a function of kp/A, normalized to the energy density £q for a free Fermi gas. The induced 
three-body contribution Eq. ( |33|) is omitted for q / 0. When it is included, all curves coincide 
with the a = curve. 



of all first order diagrams in a vanishes exactly, and the field transformation Eq. (|llj) does 
not change the energy density. 

These cancellations manifest the invariance of finite density observables under field redefi- 
nitions. It is evident that one must include all induced terms. However, in past investigations 
with unitarily equivalent two-body potentials, possible induced three-body contributions 
were not considered. (As demonstrated in the next section, unitary transformations of the 
Hamiltonian have the same effect as field redefinitions.) Fig. ^] illustrates the consequence 
of omitting these terms for a dilute Fermi gas with spin degeneracy g — 4. The energy 
density of the gas through 0(k F ) (see for details) is plotted for different values of a, 
with the contributions £e(a) an d £e(&) with two-body vertices included but the contribution 
£ 6 ( c ) with a three-body vertex not included. The result is a band of energies, analogous to 
the "Coester band" found for nuclear matter saturation curves The normalization of 
the transformation was chosen so that a = 1 is natural-sized; a manifestation is that the 
"correction" for kp = A is order unity for a = 1. The spread of the Coester line curves in 
the medium is a measure of the natural size of three-body contributions. 

In summary, we have used a simple field redefinition to illustrate how two-body off-shell 
vertices are induced but canceled exactly by new many-body vertices. Although we have 
worked to leading order in the EFT perturbation theory (that is, to leading order in a 
momentum expansion), the cancellation is guaranteed order- by-order .[] Next we show how 



As noted earlier, a similar cancellation holds in a finite system with a single particle potential 
included in the noninteracting propagator. 
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unitary transformations on NN potentials exhibit the same behavior. 



IV. UNITARY TRANSFORMATIONS AND FIELD REDEFINITIONS 

Infinitely many potentials that predict the same scattering phase shifts and energy spec- 
trum can be constructed using unitary transformations ||. If 

= (H + V)V = Em (36) 

and exp(iS) is a unitary operator of finite range: 

e iSq = $ ™ $ ; (37) 

then the transformed Hamiltonian 

H = e iS He~ iS = H + V (38) 

has the same energy spectrum and phase shifts at all energies as H. The potentials V and 
V are said to be phase-shift equivalent. Various practical methods to construct families of 
such potentials are described in the literature fl]-^|]. 

Here we demonstrate that these short-range unitary transformations have the same qual- 
itative behavior as the simple field redefinitions in the last section: off-shell contributions are 
generated at the same time as compensating many-body contributions. This is only evident 
upon embedding the problem in a many-body framework from the beginning. The interplay 
between off-shell behavior and many-body forces in the context of unitary transformations 
has been studied in the literature |]30|j31[] . In particular, Ref. [0 proves a general theorem 
that guarantees the existence of three-body potentials that, when combined with unitarily 
equivalent two-body potentials, yield identical three-body observables. 

We adapt the discussion of unitary transformations in Ref. ||31|| , which addresses the key 
points (although we will use the notation of Ref. ^0|). By writing the Hamiltonian and 
unitary operator in second-quantized form, the transformation becomes a canonical trans- 
formation of the creation and destruction operators. We consider spin-0 boson operators for 
simplicity; generalizing to spin-1/2 fermions is direct but would obscure the main points. 
We focus on a special case that highlights the relationship to local field redefinitions. 

A Hamiltonian with kinetic energy and a two-body potential V is written in terms of 
second-quantized field operators as ED| 



— *2 

H = J<Px$(x)^—$(x) + -J |rf 3 xdV^ t (x)^ t (x / )K(x,x / )^(x')^(x) . (39) 
Amghar and Desplanques |3l[] consider unitary transformations of the form 

H — ► e^He~ ld = H + i[S ,H] + ■ ■ ■ (40) 

with 

S = II rf 3 xrf 3 x> t (x)^ t (x')t/(x - x>(x')^(x) (41) 
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and U Hermitian. The simplest nontrivial form that maintains time-reversal invariance [31 
is 

C/(x) = {p,x/(x)} = i%xf(x) +x/(x)(-i)^ x (42) 

where p is the momentum operator and the curly brackets denote an anti-commutator. They 
show that the transformation (|40|)-(|42|) generates two-body off-shell contributions but also 
many-body terms. 

This transformation is nonlocal for a finite range U, which complicates the manipulations. 
If the nonlocality is of the order of the short- distance scale A, from an EFT perspective it 
is appropriate to replace U(x) with contact terms. That is, U(x) becomes proportional to 
a delta function with derivatives. The simplest choice is 

x/(a;) = aV^ 3 (x) . (43) 

where a is an arbitrary parameter. Since the interaction is singular, we can expect to 
generate divergent contributions, which we will regulate with dimensional regularization. 

The transformation is most transparent in momentum space, with field expansions (in a 
box of volume V with periodic boundary conditions) ||40|| : 



= E ^"X , ?(x) = -^E e" ik x 4 , (44) 



k V v k 

and boson commutation relations for the operators a k and 4 : 

ho 4'] = <*kk' , [a k , a v ] = [4, 4] = . (45) 
The Hamiltonian for a zero-range potential is 

H = J2^7 4«k + ^7 E 4 1 «k 2 a k3 a k4 5 k 1 +k 2 ,k 3 +k 4 , (46) 

k ZIV1 ZV kl,k 2 ,k 3 ,k4 

with Cq given in Eq. ([17]) ■ This Hamiltonian corresponds to the Lagrangian of Eq. (16) 
with the C2 term omitted. In order to obtain the transformation operator S, we substitute 
Eqs. (|42"D and (j43[) into Eq. (p|) . After changing the integrations over x and x' into integra- 
tions over z± = x±x', the z + integration simply gives overall momentum conservation. The 
remaining integral can be carried out using a partial integration, and the transformation 
operator becomes 

S = i -jb E (( k i - k 2) 2 - ( k 3 - k 4) 2 ) 4 1 4 2 ak 3 a k4 5 k 1 +k 2 ,k 3 +k4 ■ (47) 

ki,k 2 ,k 3 ,k 4 

The transformed Hamiltonian is 
H' = H + i[S,H] + 0{a 2 ) 

= ^ + 16^V ^ ((k x - k 2 ) 2 - (k 3 - k 4 ) 2 ) a kl 4 2 ak 3 a k4 5 k 1+ k 2 ,k 3 +k4 

ki,k 2 ,k 3 ,k4 



E (( k l " k 2) 2 " ( 2k 5 - k 3 - k 4) 2 ) 4 1 4 2 a k3 a k4 \ 1+ k 2 ,k 3+ k4 
ki,k 2 ,k3,k4,k 5 

M( (( k l " k 2) 2 - ( k 4 - k 5 - k 6 + k 3 ) 2 + (k 4 - V 2 



2V 2 

ki,k 2 ,k3,k4,k 5 ,k 6 



-(kx - k 2 - k 3 + k 6 ) 2 ) 4 1 4 2 4 3 a k4 a k 5 a k 6 5 k 1+ k 2 +k3,k4+k 5 +k 6 + °( a2 ) (4 
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The two-body term generated from the commutator with the potential contains a divergent 
sum over the momentum k 5 . However, this is a pure power divergence, which is zero in 
dimensional regularization (with a cut-off regulator, it would simply renormalize the Co 
coefficient in the effective Lagrangian). 

We observe that the momentum dependence of the first two-body term, which results 
from the commutator with the kinetic term, can be written 

(k! - k 2 ) 2 - (k 3 - k 4 ) 2 = 2(k 2 + k 2 - kj - k 2 4 ) . (49) 

This factor vanishes for on-shell legs and so does not contribute to two-body scattering. The 
analysis proceeds in parallel with the analysis in the last section; we find the analog of Fig. ||| 
with the same result: the off-shell and three-body contributions cancel for 3 — > 3 scattering 
or in the medium. 

We can identify a local field redefinition equivalent to this unitary transformation by 
directly transforming the field operator: 

$(x) — ► e <S V>(x)e -i ^ = #(x) + i[§, + 0(a 2 ) . (50) 



The commutator is conveniently evaluated by expanding in the momentum basis, Eq. (§4|), 
and using the commutation relations, Eq. (551) . The end result can then be interpreted as a 
field redefinition in the effective Lagrangian: 

i) — > tfj + a V)^ + ^ (^) 2 + (^ f ) ■ VW)] + 0(a 2 ) , (51) 

with the corresponding redefinition of ip\ This transformation generates a vertex equivalent 
to the term proportional to a/(MV) in Eq. (|48| ) but also appears to generate a vertex with 
time-derivatives from the term, 



AC = a 



^VWV) + (^ f )Vl (W) + H.c. , (52) 



which does not appear in the transformed Hamiltonian. However, the Feynman rule for this 
vertex is proportional to uj-l-\- uji — uj% — oj^ , which vanishes identically by energy conservation. 

In Ref. fliTf , EFT at finite density was explored by studying a perturbative matching 
example. In particular, a model potential 

(k |y true |k) = ^ (A .2 +m 2 )(Jfc ,2 + m 2) • ( 53 ) 

was used to represent the underlying dynamics, where the mass m corresponds to the range 
(and non-locality) of the potential; it plays the role of the underlying short- distance scale A. 
The dimensionless coupling 7 provides a perturbative expansion parameter. An EFT was 
introduced to reproduce observables in a double expansion: order-by-order in 7 as well as 
the EFT momentum expansion. 

An effective s-wave potential in the form 

(k'|l> EFT |k> = C + C 2 ^-±^ + C ^^P 2 + C ^-p 2 + . . . , (54) 
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was used. The EFT expansion of this potential was studied using both dimensional regu- 
larization and cut-off regularization. The term proportional to C4 does not contribute to 
on-shell two-body scattering in dimensional regularization, but is essential in the matching 
of the energy per particle in an infinite system. The value of C4 was determined by matching 



on-shell 3^3 scattering in Ref. [O]]. But it is tempting to conclude that this is merely 
equivalent to matching the T-matrices generated by Eqs. (|53|) and ( |54D off shell. 

However, the unitary transformation discussed above shows that we can arbitrarily trade- 
off the value of C4 with the coefficient of a three-body contact term. This means that off-shell 
matrix elements of the potential (e.g., (0|VEFT|k)) and subsequently the off-shell T-matrix 
can be changed continuously by varying a, without changing any observables. One choice 
is to eliminate two-body off-shell vertices such as C4 entirely in favor of many-body vertices 
that do not vanish on shell. This is the form of many-body EFT used in Ref. ||33|| , which 
corresponds to Georgi's "on-shell effective field theory" |18|| . In different situations, different 
choices may be more efficient 



V. SUMMARY 

The attitude in the mid-seventies about off-shell physics is summarized in the review 
article by Srivastava and Sprung ||: 

"A knowledge of the off-energy-shell behavior of the nucleon-nucleon interaction 
is basic to nuclear physics. The on-energy-shell information, however complete 
(which it is not), is not adequate to permit unambiguous calculation of the 
properties of systems of more than two nucleons. Any reasonable agreement with 
nuclear binding energy and other properties obtained by using different potential 
models without knowing the correct off-energy-shell behavior of the interaction 
is probably fortuitous and therefore not very meaningful. The results obtained 
with different potentials only rarely agree with each other." 

In contrast, effective field theories are determined completely by on-energy-shell information, 
up to a well-defined truncation error. In writing down the most general Lagrangian consistent 
with the symmetries of the underlying theory, many-body forces arise naturally. Even though 
they are usually suppressed at low energies, they enter at some order in the EFT expansion. 
These many-body forces have to be determined from many-body data. The key point is, 
however, that no off-energy-shell information is needed or experimentally accessible. 

In principle, one could consider the formal problem of particles interacting via a definite 
two-body potential, valid for all energies. From a physical point of view, such a description 
has to break down at a certain energy, when, for example, the substructure of the particles 
can be resolved. Furthermore, such a scenario cannot correspond to QCD where three-body 
forces are present at a fundamental level. The EFT approach is ideally suited to deal with 
such a situation, because at low energies all sensitivity to the high-energy dynamics can be 
captured in a few low-energy constants. 

The variation in observables calculated using the two-body parts only when short-range 
unitary transformations are performed gives an indication of the natural size of three-body 
contributions. Given the correspondence of unitary transformations and field redefinitions 
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(as shown in the previous section), this statement is a simple consequence of naive di- 
mensional analysis in the EFT. In cases where naive dimensional analysis fails, further 
investigation is needed. 

The modern consensus is that at least three-body forces are required as supplements 
to the best phenomenological NN interactions; discrepancies between calculated and exper- 
imental binding energies for the triton have made this conclusion unavoidable. But the 
emphasis on the importance of reproducing the "correct" off-shell physics still persists. For 
example, in a recent review of the state-of-the-art Bonn nucleon-nucleon potential, differ- 
ences in the off-shell behavior of the Bonn potential compared to other contemporary NN 
potentials are stressed |[43|| . The hope of using many-nucleon systems to decide on the "cor- 
rect" off-shell potential also persists |4^j45|| . Again, the EFT viewpoint is that all theories 



with consistent power counting must agree up to truncation errors. 

We have noted several consequences of the freedom to perform field redefinitions: 

• An effective Lagrangian with energy-independent effective potentials can always be 
used, by eliminating time derivatives with field redefinitions. However, it is sometimes 
numerically convenient to trade momentum dependence for energy dependence [|7|], or 
to retain time dependence when incorporating relativistic corrections [^2| . 



• The Lagrangian can be restricted to on-shell vertices only (cf. Georgi's "on-shell effec- 
tive field theory" [0). One advantage of this choice is that the maximal number of 
energy diagrams vanish (at zero temperature). 

• One can choose to trade off-shell contributions for many-body forces. If minimizing 
many-body forces is desirable, then one could argue that adjusting off-shell behavior 
to achieve this end is optimal. However, fine tuning is required to reduce many-body 
contributions below their "natural" size. 

Finally, we note that different regularizations will also yield very different potentials, with 
different off-shell behavior, but observables will agree. 
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APPENDIX: NUMBER OPERATOR 



The field transformation, Eq. (18), changes the number operator of the theory, which 
is defined as the time component of the Noether current associated with the phase trans- 
formation ip — > exp(—i(f))if). In the untransformed theory, the number operator is simply 
N = i[)*if). The number density no for the free Fermi gas at T = is given by the Feynman 
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(a) (b) (c) (d) 

FIG. 9. Finite Density: diagrams for the number operator. 

diagram shown in Fig. [|(a), where the cross indicates the insertion of the number operator. 
Evaluation of the Feynman diagram gives 

n = -g lim [ -^- iGo (p)e^ = g^- . (Al) 



71- 



0+7 (2tt) 4 uwv ^6tt 2 



This density n Q remains unchanged by interactions and n = n , which can be understood 
as follows. All contributions from interactions have the form shown in Fig. |](b), where the 
blob represents the remainder of a given Feynman diagram. The crucial point is that there 
are always two different propagators with the same momentum connected to the insertion 
of N. As a consequence, there is no simple pole in the time component of the corresponding 
momentum, and the contour integral over the time component vanishes. This argument is 
valid as long as there are no off-shell terms in the vertices that cancel one of the propagators 
in question. 

After the field transformation, Eq. (|18D , the number operator acquires an additional term 
of 0(a) and becomes 

^ = ^ + ^T 2 (^V) 2 • (A2) 

The expectation value of the number operator and the density, however, remain unchanged. 
We show this explicitly for the lowest order in a. At order a, we have only the two new 
Feynman diagrams shown in Figs. |9](c) and |](d). The cross indicates an insertion of the part 
of the number operator that is 0(1), while the tetragram denotes the part of the number 
operator that is 0(a). The triangle denotes the induced off-shell vertex (hi) (cf. Eq. (|20D). 
The diagram ^(c) does not vanish from the general argument above because the off-shell 
vertex cancels one of the propagators connected to the insertion of N. Evaluating the two 
diagrams, we find 

"9(c) = -9(9 - 1) A 3 36?r 4 = ' ( A3 ) 

Thus the contributions cancel and the density is unchanged to 0(a). Higher order contri- 
butions in a cancel in a similar way. 
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